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ABSTRACT 

Common models of blazars and gamma-ray bursts assume that the plasma underlying the ob- 
served phenomenology is magnetized to some extent. Within this context, radiative signatures 
of dissipation of kinetic and conversion of magnetic energy in internal shocks of relativistic 
magnetized outflows are studied. We model internal shocks as being caused by collisions of 
homogeneous plasma shells. We compute the flow state after the shell interaction by solving 
Riemann problems at the contact surface between the colliding shells, and then compute the 
emission from the resulting shocks. Under the assumption of a constant flow luminosity we 
find that there is a clear difference between the models where both shells are weakly magne- 
tized (cr< 10"^) and those where, at least, one shell has a cr> 10"^. We obtain that the radiative 
efficiency is largest for models in which, regardless of the ordering, one shell is weakly and 
the other strongly magnetized. Substantial differences between weakly and strongly magne- 
tized shell collisions are observed in the inverse-Compton part of the spectrum, as well as in 
the optical. X-ray and IGeV light curves. We propose a way to distinguish observationally 
between weakly magnetized from magnetized internal shocks by comparing the maximum 
frequency of the inverse-Compton and synchrotron part of the spectrum to the ratio of the 
inverse-Compton and synchrotron fluence. Finally, our results suggest that LBL blazars may 
correspond to barely magnetized flows, while HBL blazars could correspond to moderately 
magnetized ones. Indeed, by comparing with actual blazar observations we conclude that the 
magnetization of typical blazars is cr < 0.01 for the internal shock model to be valid in these 
sources. 
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1 INTRODUCTION 

Highly variable radiation flux has been observed in the rela- 
tivistic outflows of blazars and gamma-ray bursts (GRBs). Even 
though the radiation energy and time scales are difi"erent for both 
classes of objects (y-rays on a millisecond timescale for GRBs 
versus X-rays on a timescale of hours for blazars) the underly- 
ing physics responsible for the ener gy dissipation might be very 
similar. The internal shock scenario (IRees & Meszaroslll994) has 



been used to explain the variability of bl azars (e.g., Spada et al] 
200lUMimica et alJI 20041) and GRBs (e.g..lKobayashi et alJll997F 



Daigne & Mochkovitchlll998l : lBosniak et al.j2009l) . In this scenario 
inhomogeneities in a relativistic outflow cause parts of the fluid 
to collide and produce shocks waves which dissipate energy. The 
shell collisions are often idealized as collisions of dense shells. In 
recent years on e- and two-dimensional relativistic hydrodynamics 
(RHD. lKino"et al. 2004; Mimica et al. 2004, 2005) and relativistic 
magneto-hydrodynamics (RMHD. lMimica et alj2007l) simulations 
of the shell collisions have been performed and have showed that 
the dynamics of shell interaction is much more complex than what 
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is commonl y assumed when modeling shell interactions analyt- 
ically (e.g. Kobavashi et al.lll997l; IPaigne & Mochkovitchlll998l; 



ISpada et al.l(200ll : lBosniak et al.ll2009ir Particularly, the influence 
of the magnetic fiel d (if present) has b een shown to significantly 
alter the dynamics jMimica et alj|200"7h . In spite of these efforts, 
we still do not known with certainty whether the flow, whose en- 
ergy is being dissipated, is significantly magnetized, or whether it 
is only the kinetic energy which ultimately powers the emission. 

In a previous work jMimica & Aloll2010l . MAIO in the fol- 
lowing) we have studied the dynamic efficiency, i.e. the efficiency 
of conversion of kinetic to thermal and/or magnetic energy in inter- 
nal shocks. We found that the dynamic efficiency is actually higher 
if the shells are moderately magnetized (cr x 0.1, see the next sec- 
tion for the definition of cr) than if both are unmagnetized. How- 
ever, we did not compute the radiative efficiency of such interac- 
tions, but inst ead used the dynamic effic iency as an upper bound 
of it. R ecently Bottcher & Dermer ('2 010l BDI O in the following), 
lloshi & Bottcher (2011) and Chen et al] | |201 ih have presented so- 
phisticated models for the detailed computation of the emission 
from internal shocks. While these models assume a simple hydro- 
dynamic evolution, they employ a time-dependent radiative trans- 
fer scheme which involves the synchrotron and synchrotron self- 
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Compton (SSC) processes as well as the contribution Comptonized 
external radiation (external inverse Compton - ETC), all the while 
taking into account the radiative losses of the emitting non-thermal 
particles. We have adapted the method of BDIO and use it to per- 
form a parametric study, addressed to infer the magnetization of the 
flow from the light curves and spectra of internal shocks in magne- 
tized plasma. 

The organization of this paper is as follows: Section|2]briefly 
summarizes the model of MA 10 which is used to study the shell 
collision dynamics, and in Sections[3]and|4]we describe the numer- 
ical method we employ to compute the non-thermal radiation. We 
discuss the radiative efficiency in Section[5]and present the spectra 
and light curves in Section (6] A global parameter study is elab- 
orated in Section [7] We close the paper with a discussion of our 
results and give our conclusions (Sec. [8}. 



2 SHELL COLLISION DYNAMICS 

As was discussed in detail in the Section 2 of MAIO, our aim is to 
model a large number of shell collisions with varying properties. 
Therefore, we employ a simplified model for a single shell colli- 
sion, based on the exact solution of the Riemann problem. When 
describing the initial states of the Riemann problem we will use 
subscripts L and R to denote left (faster) and right (slower) shells, 

respectively. 

W e assume a cylindrical outflow with a radius i^. lMimica et al.l 
( I2OO4I) show that the jet lateral expansion in this case is negligible. 
For simplici ty, and being consistent with previous work in the field 
(e.g., BDlO. ljoshi & BottcheJiioTTh , we also ignore the shell lon- 
gitudinal expansion after the shocks cross the shell (see also Sec- 
tion [3]4j. Following the equation 9 of MAIO we define the lumi- 
nosity as 



(1) 



L := nR-pc' [r^(l + e+x + cr)-T\ Vl -T-^ , 

where c is the speed of light in vacuum, p is the fluid rest-mass den- 
sity, e is the specific internal energy, := p/(pc^) is the initial ratio 
between the thermal pressure and the rest-mass energy density, and 
cr ;= B^/(4npT-c^) is the magnetization parameter. Here B is the 
strength of the large-scale magnetic field, which is perpendicular 
to the direction of propagation of the fluid moving with velocity v 
and a corresponding Lorentz factor F := 1/ -sjl - (v/c)^. The spe- 
cific internal energy is related to the pressure and to the density 
through the equation of state. We use the TM analytic approxima- 
tion to the Synge equ ation of state dde Berredo-Peixoto et al .I2OO5I : 
et al]2005h and obtain: 



the particulars of the initial states these waves can either be shocks 
or rarefactions. We label the left-going wave with RS to denote a 
reverse shock, and with RR in case a reverse rarefaction happens. 
Similarly, we label the right-going wave with FS or FR to differ- 
entiate the cases in which a forward shock or a forward rarefaction 
occurs, respectively. We will use a subscript 5 to refer to the prop- 
erties of the shocked fluid in general, and the subscripts FS and 
RS when distinguishing between the front and reverse shocked flu- 
ids. Finally, we will use the subscript for properties of the initial 
states in general, and the subscripts L and R when we need to dis- 
tinguish between left and right initial states. Because we assume 
that the flow luminosity is the same for both initial states, using Q} 
we determine the number density in the shells to be 



(3) 



where Ti = (I + Ag)FR. 

The Riemann solver provides us with the bulk velocity of the 
shocked fluid/3c (and its Lorentz factor F = (1 -yS^)"''^), and veloc- 
ities ySfi-c and Prsc of the FS and RS, respectively (provided they 
exist). The velocity of the initial (unshocked) states in the CD rest 
frame is 



The shock velocities in the frame of the CD can be computed as 



(4) 



1 -ySA- 



(5) 



where prime denotes quantities in the CD rest frame. In this frame 
the shock crosses the shell at a time 



'cross.S 



where Ar'n is the shell width in the CD frame. 



A^o = FAr 



(6) 



(7) 
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3 NON-THERMAL PARTICLES 

In this section we show the properties of non-thermal particles and 
their emission. We first discuss the model for the magnetic field and 
non-thermal particles, and then outline the method used to compute 
their emission. 
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We assume that Li = Lf, and;^^^ = xr- Furthermore, as in MAIO, 
we assume '■= (1 + Ag)F/i. This leaves us with R, ctl, and 
Ag as parameters, because all other quantities can be determined 
using equations (TJ and To these, we add an additional param- 
eter Ar, which is the initial width of the shells in the LAB-frame. 
While it does not influence the solution of the Riemann problem, 
it provides the physical scale necessary for the calculation of the 
observed emission. 

Once the initial states are constructed, we com pute the exact 
solutio n of the Riemann problem using the solver of iRomero et al.l 
( l2005h . The initial discontinuity between left and right states de- 
composes into a contact discontinuity (CD), and a left-going and a 
right wave (in the frame in which the CD is at rest). Depending on 



3.1 Magnetic field 

As in iMimica & Aiovl ( I2OI0I) and BDIO, we assume that there ex- 
ists a stochastic magnetic field, which is created in situ by the 
shocks arising in the collision of the shells. We label this field by 
6s.st> and by definition its strength is a fraction eg of the internal 
energy density of the shocked shell 1*5 (obtained, in our case, by 
the exact Riemann solver): 



(8) 



Since we allow for arbitrarily magnetized shells, there is also 
an ordered (macroscopic) magnetic field component Ss.mao which 
is a direct output of the exact Riemann solver. The total magnetic 

field is then Bs 
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3.2 Injection spectrum of non-thermal particles 

We assume that a fraction of electrons in the unshocked shell is 
accelerated to high energies at the shock front. Following Sec- 
tion 3 of BDIO, we assume that a fraction of the dissipated ki- 
netic energy is used to accelerate electrons. We assume that some 
particle acceleration mechanism operates at shocks, such that a 
fraction of the electrons in the unshocked shell is accelerated to 
high energies in the vicinity of the shock front. As it is com- 
monly done I Bykov & MeszarosI 1 19961 ; lOaigne & MochkovitchI 
Il998l :l Mimica et al. 2004i, BDIO), we assume that a fraction Se of 
the dissipated kinetic energy is used to accelerate electrons. The 
width of the acceleration zone Ar^^.^. is parametrized as a multiple 
Aacc of the proton Larmor radius in the shocked fluid. Therefore, 
the acceleration at a given point in the shocked fluid lasts for a time 
A4c=A4c/C8sc).Wehave 



Ar' = A,c 



' eBs 



(9) 



From this expression, we compute the volume where acceleration 
takes place as V^^^ = nR}/S.r[^^. The energy injection rate into the 
acceleration region is 



d£' 



At' 



Ar'., 



(10) 



and we assume that the energy spectrum of the injected relativistic 
particles is a power-law in the electron Lorentz factor y. 



dn' . 
df dy 



(11) 



where n[^^ is the number density of the injected electrons, is a 
normalization factor and y\ and 72 are the lower and upper injection 
cutoffs (computed below), all measured in the shocked fluid rest 
frame. The step function is defined as usual by H{x\ a, fo) = 1 if 
a < X < and otherwise. 

A cautionary note should be added here regarding the fact that 
we choose that the spectral energy distribution of the injected parti- 
cles is a pure po wer-law even in the high -cr regime. Both theoretical 
arguments (e.g.. lKirk & Heavenslll989|) and recen t particle-in-cell 
(PIC) simulations (e.g.. lSironi & Spitkovskvl2009l) have shown that 
particle acceleration is not very efficient in the presence of a strong 
magnetic field parallel to the shock front. The modifications to the 
particle injection spectrum might involve t he presence of the ther- 
mal po pulation. Recently, a calculation by[G iannios & Spitkovskvl 
( |2009|) shows how the spectrum of the GRB prompt emission might 
look in such a case; a bump at the spectral maximum and a lower 
contribution at ultra-high energies. However, current PIC simula- 
tions have not been run for sufficiently long time to achieve a stable 
situation. Thus, the fracti on of the energy which goes in to thermal 
electrons (parameter 5 of iGiannios & Spitkovs kvl l2009l) is still to 
be determined. In this particular study we set this fraction to zero 
and thus avoid introducing another free parameter in our models. 

Integrating the distribution Eq.[TT]in Lorentz factor and equat- 
ing the result to Eq.llOldivided by V.^^^ (in order to obtain the energy 
density injection rate into the acceleration region) we can compute 
the normalization factor for the electron injection. 



dE' - Jdt' 
2o = -rrr — T X 
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\iq = 2 



(12) 



3.3 Particle injection cut-offs 

As was done in iMimica et all ( l2010h . we assume that the upper 
cutoff for the electron injection is obtained by assuming that the 
acceleration time scale is proportional to the gyration time scale. 
Then the maximum Lorentz factor is obtained by equating this time 
scale to the cooling time scale, 



72 



1/2 



(13) 



where Oacc > 1 is the acceleration efficiency parameter (BDIO). 
The lower cut-off is obtained by assuming, in complete analogy to 
Eq. [To] that the number of accelerated electrons is related to the 
number of electrons passing through the shock front, 

dA'.' 



nj.O 



dV 



(14) 



(15) 



where is the fraction of electrons accelerated into the power-law 
distribution. From Eqs.[T4l[T0]and[TT]we get 

f" dy y-'i T'gnom.c^ 

Since we are dealing with potentially highly magnetized fluids, the 
condition y2 » yi cannot be assumed (see Eq. I13t . and therefore 
we cannot use the equation such as Eq. 13 of BDIO. Therefore we 
compute yi from Eq . 1 1 5 Inumerically using an iterative procedure. 



3.4 Evolution of the particle distribution 

In this work we assume that particles cool via synchrotron and 
external-Compton processes. We ignore the adiabatic cooling in 
this work since we are primarily interested in collisions of magne- 
tized shells, where the electrons are fast-cooling. The consequence 
of not accounting for the adiabatic loses of the particle distribu- 
tion is that our method overestimates the emission after the shocks 
cross the shells. Nevertheless, most of the features that make sub- 
stantive differences between the dynamics triggered by magnetized 
and non-magnetized shells happen in the early light curve and, thus, 
neglecting the expansion of the shells plasma does not change the 
qualitative conclusions of this paper 

The radiative losses for an electron with a Lorentz factor y can 
be written as 



7- 



4 u' 

--ca-T— 
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(16) 



where u'g = /SttF^ and k^^, are the energy density of the magnetic 
field and the external radiation field (see Sec. I4.2t in the shocked 
fluid frame, respectively. Once the energ y losses have been sp ec- 
ified, we use the semi-analytic solver of iMimica et al.l ( |2004|) to 
compute the particle distribution at any time after the start of the in- 
jection at shock. More precisely, we use the so lution for the contin- 
uous injection and particle cooling (Eq. 19 of IMimica et al.l l2 004h 
for the time Af^^^.^ since the beginning of the shock acceleration. 
After that time the shock acceleration ends and we approximate 
the resulting particle distribution by a piecewise pow er-law func- 
tion. Then we employ Eq. 17 of IMimica et a i] l l2004h on each of 
the power-law segments to compute the subsequent evolution. 



4 NON-THERMAL RADL4TION 

We assume that the observer's line of sight makes an angle 9 with 
the jet axis, which is also the direction of propagation of the fluid. 
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We denote by x and t the position and time in the observer frame, 
and by x' and f the location and time in the CD frame. We assume 
that the CD is located at x' = for all f . Let yu be /j := cos 0, and 
define the time at which an observer sees the radiation emitted from 
X at time t as 



T = t-xfi/c 



(17) 



Then the time as a function of the time of observation and position 
can be written as 



f = D(T/a+z) + r(j^-/3)x'/c) , 



(18) 



where 2) := [r(l - /Jyu)] ' is the Doppler factor and z is the red- 
shift. Lorentz transformations have been applied to Eq.[T7]to obtain 
Eq.M 

An important quantity is the time elapsed after the shock has 
passed a given x'. From such value one can calculate the age of the 
electron distribution function at that position, which turns to be the 
time since the shock acceleration has begun. Thus, the age can be 
defined as 



X' 



(19) 



A more useful expression involves T. Using Eqs.[T9]and[T8]we get, 
for the FS 



X' 1 -j8f5yU 



l+z 
and for the RS 
T 



D 



1 -ySssy" 



l+Z cT liRs-P 



(20) 



(21) 



We note that Eq. |20]has to be used when x' > 0, while Eq.l21lis 
valid when x' < 0. If < 0, then the shock has not crossed that 
position yet and, consequently, that place does not contribute to the 
emission yet. 

The observed luminosity in the CD rest frame is 



v'L[,(T) = nR- J dx' v'l,[t'„(T,x')]., 



(22) 



where the lower and upper limits depend on (1) whether the shock 
exists, and (2) whether it has crossed the shelQ. If the RS does not 
exist, then x' . = 0; otherwise it is 



<„„(^) = max 



TcT Prs -P 

1 +Z I -ySaSyu' 



(23) 



where Ar[ is the faster shell width in the CD frame. Analogously, 
for x' = the FS is non-existent; otherwise it is 



™ ■ .TcT' Pfs-P . , 
AT) = mm , Ar,, 



(24) 



Ar-'f, being the slower shell width in the CD frame. We point out that, 
to perform the integral in Eq. |22l /^,(?',) should be computed for 
the particle distribution evolved using values for the FS if x' > 0, 



' Since the assumed geometry in this model is cylindrical, the effects of the 
high-latitude emission are ignored. This means that all the peaks and breaks 
in the light curves are sharper than would be in case a conical geometry 
was assumed. Another consequence of the assumed geometry is that we 
overestimate the rate at which the light cui-ve declines after shocks cross the 
shells (see Section l6^ . 



and RS if x' < (i.e. in Eqs.[T2l[T3l[T5]and[8]the values for the 
corresponding shocked fluid should be used). 

Considering that v' = y(l + z)/D, and using Eg. 1221 we can 
comp ute the flux in the observer frame obtaining (BDIO, iDermeil 
l2008t) 



yFy(T) = j dx' v'l, [t'^iT, x')] , 



(25) 



AT) 



where dt is the luminosity distance. We perform the integration in 
Eq . [25] numerically. 

The total emissivity is assumed to be the result of combin- 
ing three emission processes: (I) synchrotron radiation, (2) inverse 
Compton with an external radiation field (EIC), and the synchrotron 
self-Compton (SSC) up-scattering. These emission processes are 
considered in more detail in the next sections. 



4.1 Synchrotron emission 

We compute the synchrotron emission for each power-law segment 
of the electron distribution (see Sec. I3.4t separately. In order to 
sp eed up the calculation we use the interpolation method described 
in iMimica et alj ( |2009| . section 4) and, in more detail, in I Mimical 
( |2004 sections 2.1.3 and 4.3.1). 



4.2 EIC emission 

Following BDIO, we assume that the external radiation field is 
monochromatic and isotropic in the AGN frame. We denote the fre- 
quency and the radiation field energy density in this frame by Vext 
and i<ext, respectively. Transforming into the shocked fluid frame 
we get 



FVext 
F^K,, 



(26) 



Analogously to the computation of the synchrotron emission 
(Sec. 14. It . we compute the EIC emissi vity for each po wer-law seg- 
ment separately. We use Eq. 2.94 of iMimic 1 ( |2004|) . but replac- 
ing /(vo)/vo by cM^x,/v^x, and with an additional cut-off (approx- 
imating the Klein-Nishina decline of the Compton cross-section) 
such that the emissiv ity is zero for hv > m^c* KhV^^^ (see also, 
IaIov & Mimical2008l) . Values of Vext and «ext used in this work can 
be found in Table 1. 



4.3 SSC emission 

Analogously to Sec.|4!2l we use the Eq. 2.94 of [Mimical ( |2004|) . 
However, in the case of SSC the incoming intensity of the syn- 
chrotron radiation depends on x' and T. For a point on the shell 
axis the (angle averaged) intensity at frequency vo can be written as 



7: I 

/o,v„(r,x') = ^ J d0' J d.;-;. 



„(f'(r) - sjc, x' + s cos (?') , (27) 



where L{6') is the length of the segment in direction 9' from which 
synchrotron emission has had time to arrive to x' at a time t' , and 
t'(T) is computed using Eq.[T8] The synchrotron emissivity 
can be rewritten in terms of T using Eqs.[T9] andl20l(orl2It. 



fJ.T,x')- 



1 + cos 9' 



(28) 
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Parameter 


value 


Tr 


10 


Ag 


1 


0"L 


[10-6,10'] 


O-R 


[10-6, igl] 


€B 


10-3 


€e 


10-1 


ie 


10-2 


AgcC 


10 


flacc 


10* 


R 


3 X 10'6 cm 


Ar 


6x 10'^ cm 


9 


2.6 


L 


5 X 10"^ erg s- 


"ext 


5 X 10-* erg cm 


Vext 


10''' Hz 


Z 


0.5 


e 


5° 



Table 1. Blazar model parameters used in this work. Note that at and (Tr 
can vary continuously in the indicated range. 

From Eq.[28]we can see tliat L{6') can be computed by requiring 
tliat the following condition be satisfied for each ff , 

f;(r,x')--fi+cose'i-:^l>o. 

If this condition is not satisfied, it means that the shoclc has not 
passed the point x' + s cos 9' at time t'(T)- s/c yet , i.e. there is no 
synchrotron emission from that point to contribute to the incom- 
ing intensity. In addition, we also require that L(9') < R. Finally, it 
should not be forgotten that when x' + s cos 6' > the emissivity 
of fluid shocked by the FS should be used, and the one correspond- 
ing to the shoclied fluid by the RS otherwise. Also, if either of the 
shocks is not present, there is no contribution from the correspond- 
ing region. 

In practice the numerical cost of a direct evaluation of double 
integral in Eq. [27] is prohibitive if we take into account that this 
intensity has to be evaluated for each x' in Eq. [25] To overcome 
this problem we approximate Eq. [27] by discretizing the angular 
integral in a non-uniform 6' -intervals. The choice of non-uniform 
intervals is motivated by the fact that most of the contribution of the 
incoming radiation comes from angles close to yu = -/3s , so that we 
concentrate most of the bins close to that angle. Numerical testing 
shows that using 13 bins provides an acceptable tradeoff between 
the accuracy and the computational requirements. 



5 RADIATIVE EFFICIENCY 

In this section we compare the radiative efficiency of the inter- 
nal shocks with their corresponding dynamic efficiency. We use 
the Idnematic parameters from MA 10 in the blazar regime, while 
the parameters used to compute emission are guided by the values 
from BMIO (see Table [T] for the complete list). All parameters of 
our models are fixed except for <Ti and cr^, which can vary in the 
range indicated by the Table [T] In the rest of the paper we distin- 
guish models by the value of the magnetization of each shell, e.g., 
a model with cr^ =0.1 and a-jf = I is denoted by the pair (0.1, 1). 

As can be seen from Eq. \T0\ in our model only the thermal 
energy can be injected into the non-thermal particle population. 



We point out that, alternatively or simultaneously, magnetic dis- 
sipati on can provide a source for the emission in internal shocks 
(e.g.. iGiannios et al]|2009l : iNalewaiko et al.ll201 ll . and references 
therein), which we are not considering here. Thus, the radiative effi- 
ciency we compute in this paper is only a lower bound to the actual 
radiative efficiency of the binary collision of relativistic magnetized 
shells. As is shown in the Appendix IaI we can use t he definition of 
the dy namic efficiency inspired by the recent work of lNaravan et al.l 
( |20T ij). Its advantage is the Lorentz-invariance, which enables us to 
compare it to the radiative efficiency of our model. 

Following MA 10, but using the definitions for the different 
energy components E^, Et and Em of the Appendix (Eq. lA2b . we 
denote by Eq the total energy in the unshocked shells, 

Eq := Ek(Tr(1 + Ag), ntm,,, 1) + Et(Tr{1 + Ag), nLmp,xnLmpC^ , 
+Em(Xr{'^ + ^g),iiLmp,(rL, 1) + EK(TR,iiRmp, 1) 

+Et{T R,nRmi„xnRm,,c^ , 1) + EmiTR^nRmp^aR, 1) . 

where is the pressure-to-density ratio of the cold initial shells and 
is set to 10 We also define the width of the shocked shells in 
terms of their initial width, 

r .-\izM no^ 

so that ^ < 1 H. The dynamic thermal efficiency for the faster 
shell is defined as 
1 

er.L := —X 

[ET(^s,L,ns,Lmp,ps,L,^w,L) - Et{Tr(\ + Ag),ni^mp,xnLmpC^, 1)J 

and analogous definitions can be written for €m.l, ^t.r and £m.r (see 
Eqs. 13, 14, 16 and 17 of MAIO). The total (Lorentz invariant) 
dynamic thermal and magnetic efficiency is 

gj. = gjj^ + , (32) 

fM = fM.L + £m,« ■ (33) 

From these equations it can be seen that the radiative efficiency can 
be at most eder + Sm)- More formally, we can write the radiative 
efficiency as (neglecting adiabatic cooling) 

fiad := ee/rad(er + <?«) , (34) 

where /^ad := er/ier + ^m)- It should be noted that Eg. 1341 refers 
to the "bolometric" emission, i.e. it includes all frequencies for the 
whole duration of the shell interaction. Since Earth-based observa- 
tions have a limited spectral and temporal coverage, the Eq. |34|is 
only an upper limit for radiative efficiencies inferred from actual 
observations. Figure [T] shows that the radiative efficiency is not a 
one-to-one map of the total dynamic efficiency. In particular, we 
note that /^d drops to under 10% in the region (ct/. > 10, ctr > 10). 
Furthermore, there is a region of maximal dynamic efficiency for 
ctr x 0.2 and ct;, > 1, where the radiative efficiency from inter- 
nal shocks, which can only tap the thermal energy in the regions 
downstream shocks, is not maximum. Nevertheless, for small-to- 
moderate values of the magnetizations of both shells (lower right 
quadrant of Fig.[T](, the radiative efficiency is a good proxy of the 
total dynamical efficiency. 

2 Note that comparing Egs.l7land l30l we see that Ar^ = F^n/Ar, which is 
just a Lorentz transformation of the shocked shell width from lab to CD 
frame. 
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Figure 1. Contours: /jaj (see Eq. 1341 for different values of 
{(Tl,(Tr). The contours indicate the following values of /jaj (per cent): 
1,5, 10, 20, 30, 40, 50, 60, 70, 80, 90, and 100. In the region of the param- 
eter space above the dashed line there is no FS, while the RS is always 
present for the considered parametrization (see also Fig. 1 of MAIO). Filled 
contours: total dynamic efficiency ej + £m in per cent. 



6 SPECTRA AND LIGHT CURVES OF MAGNETIZED 
INTERNAL SHOCKS IN BLAZARS 

Our aim is to produce synthetic spectra and light curves from 
our numerical models of the interaction of two relativistic, mag- 
netized shells. With this purpose, we chose three models from 
our parameter space, which are representative of different condi- 
tions that can be encountered in blazar jets. The first model cor- 
responds to a regime of very low magnetization of both shells 
(ctl, (Tji) = (10^*, 10"''). The second and third models correspond to 
intermediate (10"", 10"^) and moderate/high shell magnetizations 
(1,0.1). 



6.1 Average spectra 

The spectrum of the weakly magnetized model (10"', 10"'), (Fig.[2j 
full lines) reproduces the typical double-peaked spectrum of 
blazars. The synchrotron emission (Fig.[2j solid red line) peaks at 
4.6 X 10'^ Hz, while the inverse-Compton (IC) emission, dominated 
by the SSC component, peaks at 6.7 x 10^' Hz (Fig. [2j solid blue 
line). In this case, the IC spectral component is clearly dominat- 
ing the overall spectrum. At intermediate magnetizations (Fig. [2j 
dashed lines) the synchrotron emission peaks at higher frequen- 
cies than in the previous case, namely, at 7.5 x lO''' Hz, while the 
IC emission peaks at 6.0 x 10^" Hz, as one would expect, since 
a larg er magnetic fi eld increases the synchrotron peak frequency 
(e.g., lMimica|[2004l) . For these shell magnetizations, the SSC also 
dominates the high energy emission, but now both the SSC and EIC 
components are significantly weaker than in the model (10"', 10"'). 
Interestingly, moderate to high magnetizations (Fig. [2j dot-dot- 
dashed lines) reduce the peak frequency (1.9 x 10'* Hz) and sub- 
stantially flatten the synchrotron spectral component w.r.t. the inter- 




Figure 2. Averaged spectra (time integration interval: 0-100 ks) for mod- 
els (10"', 10"'), (10"2, 10"2) and (1,0.1) (full, dashed, and dot-dot-dashed 
lines, respectively). Black colored lines show the total spectrum, while the 
red, blue and green lines show the contribution due to synchrotron, SSC and 
EIC emission, respectively. 



mediate magnetization case. Furthermore, IC spectral components 
are notably weaker than the synchrotron one for large magnetiza- 
tion, being the IC spectral peak located at 4.2 x 10'^ Hz. 

As can be seen from Fig. [21 the synchrotron emission from all 
three models is of comparable intensity, while the IC emission is 
much weaker in the strongly magnetized model (1, 0.1). The reason 
for such a large difference in the high-energy emission between the 
magnetized and the non-magnetized models lies in the lower num- 
ber density of emitting electrons (Eq.O and in the higher magnetic 
field of the magnetized model. The magnetized model has a much 
lower number density due to its relatively high cr, in both FS and 
RS emitting regions, which means that there are less scatterers for 
the SSC and EIC processes. A high magnetic field also means a re- 
duction of the upper injection cut-off (Eq.ll3ll, which in turn means 
that the seed synchrotron photons in the SSC process are being up- 
scattered to lower frequencies, explaining the small contribution of 
the SSC component to the spectrum. The EIC component's upper 
cut-offs are determined by the Klein-Nishina decline (see Sec. |4.2| ( 
and not by the upper injection cut-off, which explains why the EIC 
spectral peaks of the models are in a similar frequency range. 

Figure[3]shows the contributions from the FS and RS (red and 
blue lines, respectively) to the total spectrum (black lines). Except 
close to the local minima between the two spectral peaks, the RS 
contribution is dominant in the average spectra of the models with 
low and intermediate magnetizations, (10"', 10"') and (10"^, 10"-) 
respectively. In the vicinity of the aforementioned minima (located 
in the X-rays range), the FS contribution tends to broaden the width 
of the minima and to soften the spectral slope. At the moderate to 
high magnetizations of the model (1,0.1) the FS is dominant except 
in the range 10'^ - 10" Hz, where the FS and the RS have com- 
parable contributions. The reason is that the faster shell, through 
which the RS propagates, is substantially more magnetized than 
the slower shell, so that the RS has less particles to accelerate than 
the FS. 
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Figure 3. Same as Fig. [2] but distinguishing the contributions of the FS (red 
lines) and of the RS (blue lines) to the total spectrum (black lines). 
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Figure 4. Multi- wavelength light curves for the models (10 ^,10 
(10^^, 10"^) and (1, 0. 1), shown with full, dashed and dot-dot-dashed hnes, 
respectively. The R-band (5 X lO'^Hz), X-ray band (1 - lOkeV), as well 
as 0. 1 GeV and 1 GeV light curves are shown in black, red, green and blue, 
respectively. 



6.2 Light curves 

The multi-wavelength light curves of the models presented in the 
previous section are displayed in Fig.|4] We have picked up several 
characteristic bands to analyse the data (R-band, X-ray, 0. 1 GeV 
and 1 GeV). 

Comparing the R-band light curves of the three models we 
can see that models (10^^, 10"^) and (1,0.1) exhibit properties of 
the fast-cooling electrons emitting synchrotron radiation, while the 
model (10^*, 10"*) shows the opposite, slow-cooling behavior. In 
the latter case the maximum of the R-band light curve is reached 
when the shocks cross the shells, and afterward the emission de- 
cays as the particles cool down (no new particles are accelerated 
after both shocks cross the shell). In the case of the model (1,0.1) 
one can clearly notice two sudden drops in emission around 4ks, 
which correspond to the moments when first the RS, and later, the 
FS cross their respective shells. The almost vertical drops in emis- 
sion are indicative of a very efficient electron synchrotron-cooling 
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Figure 5. X-ray light curves for the models (10"^ 10"'^), (10"^ 10"^) and 
(1,0.1), shown with full, dashed and dot-dashed hnes, respectively. The to- 
tal light curve is shown in black, while the synchrotron and SSC contribu- 
tions are shown in red and blue, respectively. 
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Figure 6. 1 GeV light curves for the models (10"'^, 10"*), (10"^, 10"^) and 
(1,0.1), shown with full, dashed and dot-dot-dashed lines, respectively. The 
total light curve is shown in black, while the EIC and SSC contributions are 
shown in red and blue, respectively. 



1^ At intermediate magnetizations (10"^, 10"^) the first sharp drop 
is observable as well, though here there is a weak late-time optical 
emission between 10'* and 10^ seconds due to the SSC process. 

The emission in the X-ray band is a bit more involved, and to 
perform a proper analysis we show in Fig. [5] both, the total light 
curve (black lines), and the individual contributions to it of the syn- 
chrotron and SSC processes (red and blue lines, respectively^. Ex- 
cept at very early times, the emission is dominated by the SSC pro- 
cess in all cases. The synchrotron emission in this band happens 



^ Note that, since the high-latitude emission is ignored due to cylindrical 
geometry, the drops are too sharp and would be smoother were a conical jet 
geometry assumed. 

We do not show the EIC light curve because its contribution is negligible 
at these frequencies. 
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Figure 7. Contours of the logaiithm of the time (0 - lOOks) and frequency 
(lO'^ - 10^^ Hz) integrated flux (i.e., fluence) as a function of the shell 
magnetization cr^ and cru. Note that, different from Fig. [T] the region of 
ultra-high magnetizations (<tlh > 10) is not included in this figure. The 
reason being that the computation of the integrated flux with such extreme 
magnetic fields requires a discretization of the two-dimensional integral in 
Eq.|22]in a very large number of intervals, making the calculation numeri- 
cally impractical. Nevertheless, the trends at such high magnetizations can 
be easily extrapolated from the values displayed in the figure. 



in an efficient fast-cooling regime, which can be inferred from the 
fast drop of the synchrotron components between 4 and 9 ks. The 
fact that increasing magnetic fields make that particles cool faster, 
explains that the non-magnetized model peaks much later 60 ks) 
in this band than the other two (more magnetized) models. 

At energies of 1 GeV, there is only emission from IC pro- 
cesses (Fig. [6]l. The model with the smaller magnetization dis- 
plays a clearly dominant EIC emission at early times, while in 
the other two models EIC dominates the later times. In the mod- 
els (10^-, 10"") and (1,0.1) EIC, similar to the synchrotron emis- 
sion in the X-ray band, sinks very quickly before 8 ks, indicating 
that the electrons are in a fast-cooling regime. In the latter mod- 
els, because of the delays associated to the physical length of the 
emitting region, the SSC contribution peaks very early and decays 
exponentially before the sharp drop of the EIC emission (this is par- 
ticularly the case of the most magnetized model, in which the SSC 
component does not significantly contribute to the light curve after 
^ 400 s). In contrast, the EIC emission of the model (10"*, 10"*") 
shows a much more prominent peak and a shallower decay from 
the maximum (at =^ 9 ks), both features being characteristics from 
electrons in a slow-cooling regime. 



7 GLOBAL PARAMETER STUDY 

In the following we present the results of the global parameter study 
of the dependence of the emitted radiation on the shell magnetiza- 



tior[f|. Figure|7]shows the fluence as a function of ctl and ctr. We can 
see that, as expected, the fluence roughly follows fji (see Sect.|5](. 
The region with most luminous internal shocks (upper left comer 
of Fig. |7]( happens for a moderately-to-strongly magnetized slow 
shell and a weakly magnetized fast shell, whereby the FS does not 
exist. The emission weakens as the magnetization of the fast shell 
increases, with the exception of the region where the fast shell is 
strongly magnetized but the slow shell is weakly magnetized (lower 
right corner of Fig.|7}. We conclude that, as was indicated in the 
Section 4.4 of MAIO, a large difference in the magnetization of 
the shells yields stronger dissipation and more luminous internal 
shock(s) than when both shells are weakly magnetized. 

Figure [8] shows the spectral maxima of the synchrotron, 
Vmax.syn, and the inverse Compton emission, Vmax.ic (left and right 
panels, respectively). From Fig. (2] one could anticipate a trend we 
confirm here, with the parametric space coverage, namely, that the 
IC emission is more sensitive to changes in the magnetization than 
the synchrotron emission. This statement reflects itself in Fig. [8] 
through the fact that the range of variation of v,nax,ic is larger than 
that of Vniax.syn- Thus, the IC spectral peak becomes a better proxy 
of the magnetization of the shells than the synchrotron peak. Ex- 
cept at small shell magnetizations, the IC emission happens in a 
fast-cooling regime, and the dependence of Vmax.ic with the mag- 
netization is similar to that of Complementary, at small shell 
magnetizations, the map of v„,.„.syii> resembles very much to that of 
/rad (compare the lower half of Figs.[T]and[8]-left). 

The left panel of the figure|9]shows the ratio of the IC to syn- 
chrotron fluence. The trend is quite similar to that of the integrated 
flux shown in Fig. |7] When both shells are strongly magnetized 
(cr>0. 1) the IC emission drops significantly. In the region ctl < 10"' 
the ratio is between a unity and ^ 60, with a similar behavior in 
the region (cti > 0.1, ctr < 0.1). The region of low radiative effi- 
ciency around tr^ ^ 0.01 appears as a dark vertical band in the 
plot {cTjt < lO""*, where both synchrotron and IC processes provide 
a similar fluence. The right panel of Fig. |9] shows the ratio of the 
frequencies of the IC and synchrotron spectral maxima shown in 
Fig. [8] From an observational point of view, the fluence might be 
much more robust and significant than the peak IC and synchrotron 
frequencies, which can be difficult to measure. However, for the ra- 
tio vic/vsy„, the lower right and upper left corners of the plot display 
noticeably different values. Thus, one can use the frequency ratio 
together with the fluence ratio in order to break the degeneration in 
the fluence ratio when one of the shell is very magnetized and the 
other is not magnetized (see Sec. l8.2l for further discussion of this 
point). 



8 DISCUSSION AND SUMMARY 

We have extended the study of the dynamic efficiency performed in 
MAIO by computing the multi- wavelength, time-dependent emis- 
sion from internal shocks. In this section we discuss and summarize 
our findings. 



^ All two-dimensional plots in this section have been produced using a 
logarithmically spaced grid of 30 X 30 in the o"i X rr^ parameter space. For 
each of the points we computed light curves on 96 logarithmically spaced 
frequencies (between lO'^ Hz and 3 X 10^' Hz) for 120 logarithmically 
spaced points in time (between 2 and 10^ seconds). A finer coverage of the 
<tl X a-R parameter space was not practical due to the prohibitively high 
memory and computational time requirements on the available machines. 
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Figure 8. Contours of the frequency of the spectral maxima of the synchrotron (left panel) and of the inverse Compton (right panel) emission as a function of 
o"i and ctr. 
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Figure 9. Left panel: contours of the logarithm of the ratio of the IC and synchrotron fiuence as a function of the shell magnetization cri and tr^;. Right panel: 
same as left panel, but for the ratio of the frequency of the spectral maxima of the synchrotron and inverse Compton emission. 



8.1 Emission mechanisms and magnetization 

In Section [6] we show the average spectra and multi-wavelength 
light curves of three typical models from the parameter space con- 
sidered in this paper. Synchrotron emission dominates for v < 10" 
Hz, and is rather independent of the shell magnetization (Fig.O. 
The RS dominates synchrotron emission for weakly magnetized 
shells, while in the case of strongly magnetized shells the FS and 
RS have comparable contributions. R-band light curves (black lines 
on Fig. |4]l show that the synchrotron emission is due to the slow- 



cooling electrons only for the weakly magnetized model, while for 
shells with cr > 0.01 electrons are fast-cooling. 

The SSC emission dominates in the X-ray band and higher fre- 
quencies (Fig.[2j. However, at early times the synchrotron emission 
dominates in X-rays (Fig.|5}, while in y-rays the situation is more 
complex. For the weakly magnetized model (slow-cooling elec- 
trons), EIC dominates the early emission, while in the moderate- 
to-highly magnetized models EIC dominates the late-time emis- 
sion. The reason for this is that, in the magnetized models, the 
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high-energy tail of the electron distribution disappears very quickly, 
so that the incoming synchrotron photons cannot be up-scattered 
into the 1 GeV range. In the weakly magnetized models there are 
enough slow-cooling electrons at sufficiently high energies for the 
SSC to dominate over EIC at later times. 

Finally, from Fig.|2]we see that the IC emission is weaker the 
more magnetized the shells. This is due to the requirement of our 
model that the shell luminosity (Eq.[TJ be constant regardless of a. 
From Eq.[3]we see that for cr » 1 the number density in the shells 
behaves as « cr"'. Since the IC emission depends on the number of 
electrons (EIC linearly and SSC quadratically), it is clear that the 
IC emission muss necessarily drop for large cr. From the analysis 
of the three representative models, we conclude that the shell mag- 
netization imprints two main features on the emission properties of 
blazars. On the one hand, the magnetization changes the ratio of 
integrated flux below the synchrotron peak to the integrated flux 
below the IC-dominated part of the spectrum. On the other hand, 
the magnetization determines whether electrons are slow-cooling 
(for weakly magnetized shells) or fast-cooling (moderate-to-high 
magnetization). 

8.2 Global trends 

We performed a global parameter study (Section [TJ to investigate 
the dependence of some observational quantities on the shell mag- 
netization. As discussed in Section [7] the integrated flux (Fig.|7j 
follows the trend already shown by the radiative efficiency (Fig.[Tll. 
However, the integrated flux and the spectral maxima are quanti- 
ties dependent on the particular values we have taken in our model, 
specifically, on the physical size of the shells and their bulk Lorentz 
factors, as well as on the source redshift. On the other hand, in 
MA 10 we show that the dynamic efficiency is very weakly depen- 
dent on the shell bulk Lorentz factor, i.e. it only depends on the 
shell magnetization for a fixed Ag. In order to eliminate the depen- 
dence on absolute quantities in Fig.|9]we show IC-to-synchrotron 
flux ratio, as well as the ratio of the frequency of IC to the syn- 
chrotron spectral maxima. The shape of the contours on the left 
panel of Fig.|9]does not exactly follow the one in Fig.[T] there is a 
much stronger dependence on cr^ than on cr^ in most of the scanned 
parameter space. Nonetheless, in the lower half of the plots, the be- 
havior of both Fic/Fsyn and /r is similar. For example, if we keep 
CTn constant and equal to 10"* and vary a^, we note that the ra- 
diative efficiency is larger than 90% for cr^ < 0.01, then it decays 
to a local minimum, and successively grows again to reach values 
in excess of 90% (cr^ ^1). Comparatively, at small values of cr^ 
the ratio Fic/Fsyn is close to its maximum for 0.01 ^ cr^ 1, and 
touches a minimum in the same interval as /r. The upper half of 
Figs. [9] (left) and [T] does not show the same qualitative behavior. 
The reason for this discrepancy is that RS dominates the emission, 
and thus the overall radiative properties are more sensitive to the 
magnetization of the fast shell through which it propagates. 

Interestingly, there is a certain degree of degeneration in the 
values both of the radiative efficiency and of the Fjc/Fsyn ratio 
considering the regions where one of the two shells is very mag- 
netized and the other is basically non-magnetized (i.e., the upper 
left and lower right comers of Figs.[9]and[B. In both cases, the ra- 
diative efficiency and the fluence ratio are close to their respective 
maximum values. However, we can distinguish between the case of 
high-cr/,/low-cr^ and the case of low-cr^/high-cr/i by looking at the 
ratio of peak frequencies Vic/Vsyn (right panel of Fig.|9}- A notice- 
ably smaller Vic/Vsyn ratio corresponds to the former case than to 
the later. 



The previous analysis suggests that with the combined infor- 
mation of the fluence and peak-frequency ratios, one could try to 
figure out, by using observational data, which is the rough magne- 
tization of the shells of plasma whose interaction yields flares in 
blazars. To serve such a purpose, we display in Fig.[TO]our models 
in a 2D parameter space whose horizontal and vertical directions 
are vic/Vsyn and Fic/F^y^, respectively. We notice that the computed 
models are distributed in a broad region which, nevertheless, shows 
a relatively tight correlation between Fic/F^yn and vic/Vsyn. In the 
left panel of Fig.[TOl we display our models in three different colors 
according to the magnetization of the left shell. The same has been 
done in the right panel, but for the right shell. 

Based on the degree of variation of magnetization between 
the fast and the slow shells, we have divided the parameter space in 
three broad regions (labeled with roman numerals I, II and III in 
FigllOb. where the shells have the following characteristics: 

I: moderately magnetized fast shell colliding with a weakly 
magnetized slow shell, or weakly magnetized fast shell inter- 
acting with a strongly magnetized slow one; 
II: strongly magnetized fast and moderately magnetized slow 
shells; 

III: strongly magnetized fast and slow shells. 

The first thing to note is that for models in the region III both the 
IC emission and its frequency maximum are lower compared to the 
rest of the models. This leads us to the conclusion that when the 
flow is strongly magnetized and the magnetization does not vary 
substantially the IC signature is expected to be relatively weak. Fur- 
thermore, region II shows that in the case of a larger variation in 
magnetization (i.e., weakly magnetized slow shell) the frequency 
maximum remains low, but the IC signature becomes substantially 
higher. Finally, in the region I we see that when the variation in 
magnetization is more extreme (i.e. a collision of a weakly and a 
strongly magnetized shells) we get a very strong IC signature and 
its frequency maximum is shifted to much higher energies. 

8.3 Radiative and dynamic efHciency 

As discussed in Section [5] the radiative efficiency e(,/ad(fr + Cm) 
does not have a one-to-one correspondence to the dynamic effi- 
ciency (et + fiw)- While the latter peaks in the region cr^ a 1 
and CTj; X 0.2, the former reaches its maximum in the region 
{ctl $ 10 '*,cr;i > 10). The same can be concluded from the time- 
and frequency-integrated flux shown in Fig. [7] For purposes of the 
rest of this discussion we will use /.jj as a proxy for the radiative 
efficiency. 

We note that in the region of maximum ffie FS does not 
exist. However, we see another region of high /j^j in the opposite 
corner of Fig.[T]the efficiency is quite high as well. Consistent with 
the discussion in the previous subsection and with what is shown 
in region I in Fig. [TO] we conclude that the radiative efficiency is 
maximal when the variation in magnetization between the colliding 
shells is large. 

8.4 Conclusions and future work 

Under the assumption of a constant flow luminosity we find that 
there is a clear difference between the models where both shells are 
weakly magnetized (cr < 10"^) and those where, at least, one shell 
has a cr > 10"^. We obtain that the radiative efficiency is largest in 
those models where, regardless of the ordering, there is a large vari- 
ation in the magnetization of the interacting shells. Furthermore, 
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Figure 10. Left panel: fic/fsyn (ratio of the IC to the synchrotron fluence) as a function of V[c/Vsy„ (ratio of the maximum spectral frequencies in IC and 
synchrotron ranges) for the models considered in Sec.|2] The models have been grouped in three bins according to cr/. and are annotated with black circles 
(10~' < ai < 10"*), blue triangles (10"'* < cr/. < 10"^) and red diamonds (10^^ < tr^ < lO'). Right panel: same as left panel, but in this case, the models 
have been grouped according to ctr m the same bins and same coloring and symbols in as the left panel. Shaded areas denote three regions of interest (see text 
for details). The three reference models in this paper are marked with purple crosses, and their respective magnetizations are overlaid. 



substantial differences between weakly and strongly magnetized 
shell collisions are observed in the inverse-Compton part of the 
spectrum, as well as in the optical. X-ray and 1 GeV light curves. 

In the previous sections we have deepened our analysis of the 
radiative efficiency of the process of collision of magnetized rel- 
ativistic shells of plasma. We have studied this problem from a 
mostly theoretical point of view, where the intrinsic properties of 
the flow (in particular the magnetization) have been related to the 
properties of the resulting (synthetic) spectra and light curves. It is, 
however, worthwhile to provide suitable links between our theoret- 
ical results and the observed properties of blazar flares. Thus, we 
propose a way to distinguish observationally between weakly mag- 
netized from magnetized internal shocks by comparing the max- 
imum frequency of the inverse-Compton and synchrotron part of 
the spectrum to the ratio of the inverse-Compton and synchrotron 
fluence. 

For a given flare taken in isolation, our model may predict 
which is the range of magnetizations, which have to be invoked to 
fit the spectrum. However, such a model fitting is not fully satisfac- 
tory, since it is strongly dependent on the details of the theoretical 
model. A more generic knowledge of the physical conditions in 
the flaring regions can be obtained by arranging the observational 
data in plots where the fluence ratio Fic/Fsy„ is represented versus 
vic/vsyn- The reason being that the fluence and frequency ratios are 
redshift and source independent, since they are mostly influenced 
by the variation of the bulk magnetization of the blazar jets (assum- 
ing that the viewing angle is fixed). We note that different flares of 
the same blazar, as well as different flares of distinct blazars can be 
plotted in such graphs and compared with our theoretical predic- 
tions. In addition, an average over a number of flares of the same 
source might also be interpreted using our model if the magnetiza- 
tion ratio of different pairs of colliding shells is similar. 

Our results suggest that the variability in the flow magneti- 
zation is a factor that shall be considered to explain the observed 
continuity of properties of the blazar sequence (e.g. , Fossati et al.l 
ll998l : lGhisellini et al.|[l998l : lOhisellini & TavecchidboOSi) . Look- 
ing at Figs. [2] and [3] it is evident that if the magnetization of the 



shells is not too large (at^R ^ 10"^), increasing the flow magneti- 
zation shifts the synchrotron peak towards the UV band and low- 
ers the IC peak. Leaving aside orientation effects (which we are 
neglecting here since we fix the viewing angle), such a behavior 
suggests that LBL blazars may correspond to barely magnetized 
flows, while HBL blazars could correspond to moderately magne- 
tized ones. If the magnetization is large (ctl^r ^ 0. 1 ) the synchrotron 
peak shifts towards lower frequencies, and the IC spectral peak falls 
three orders of magnitude below the synchrotron peak. The latter 
situation seems not to be observed and, thus, we conclude that this 
is an indication that the typical value of the magnetization in the 
flow of blazars is a< 10"-. We note that this value is about one order 
of magnitude smaller than that suggested for the flow in gamma-ray 
bursts (e.g., Giannios & Spruit 2006). 

Our results are not in contradiction with the common view, 
according to which, the variation of the properties of the blazar 
family correspond to the changes in the bolometric lum inosity of 
the synchrotron component (Lhni jFossati etal ]| 19971) . In such a 
case LBLs and HBLs are extrema of a one-parameter family, where 
LBLs are more radio luminous than HBLs. However, L^oi sy„(. is not 
a "direct" physical property of the plasma in a relativistic jet, since 
the same Lboi.sync may arise with an infinite number of combinations 
of bulk Lorentz factors, blazar orientations with respect to the line 
of sight, flow magnetization, etc. In this paper, we explore the role 
that the flow magnetization (a direct physical property of the emit- 
ting plasma) plays in shaping not only Lboi.synci but also the whole 
spectrum. We conclude that the flow magnetization alone is enough 
to explain the difference in Z-boi.sync and in the high-energy part of 
the spectrum (i.e., in the Compton- dominated regime) f ound in the 
blazar family. We also point out that lFossati et al.|jl998h also arrive 
to the same qualitative conclusion, since they explain that fixing the 
bulk Lorentz factor, and assuming that the SSC model be valid for 
all sources, the spectral differences in the blazar sequence shall re- 
sult from a systematic variation in magnetic field strength, HBLs 
having the highest random field intensity. We go in this paper an 
step further. When considering the dynamical changes induced by 
non-negligible macroscopic magnetic fields, the same conclusion 
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as in lFossati et alj ( fl998h holds, but we remark that, in our case, the 
total jet luminosity is kept constant (by construction of our models), 
only the magnetization is varied. 

A problem with the internal shock scenario is that it is appar- 
ently not able to explain the ultra-fast variabil ity of TeV blazars 
(e.g., [Albert et alJboOTi : | Aharonian et al. I l2007h . In order to prop- 
erly account for this f ast variability (on the time scales of min- 
utes) fast "min ijets" (Gianni os et al.l 1 20091) or "spines/needles" 
jTavecchio & Ghisellini 2008,) need to exist in a much larger, 
slower jet. We do not consider minijets in our current models, al- 
though they will need to be considered in the future. 

In the future work we will improve our modeling by includ- 
ing the resistive dissipation as a source of energy for the radiating 
non-thermal particles. This will make it possible to asess whether 
radiative and dynamical efficiencies have a one-to-one correspon- 
dence or if there is a fundamental degree of independence. In the 
latter case it might be difficult to infer the flow properties except in 
those asymptotic cases where the dissipation either does not play 
a significant role or it dominates the dynamics. Furthermore, we 
will study how changing the viewing angle and the Doppler fac- 
tor reflect on the observational signature. Finally, we will include 
the effects of the presence of a therm al component during parti- 
cle injection at very magnetiz ed shocks jSironi & Spitkovskvl2009l : 
lOiannios & SDitkovskvll2"00^ . 
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APPENDIX A: COMPARISON OF THE DYNAMIC 
EFFICIENCY DEFINITIONS 

In this appendix we compare two alternative definitions of the 
dynamic efficiency, one due to MAIO and the other inspired by 
iNaravan et al] ( 1201 iL NKTl I in the following). 

MAIO define in their Eq. 12 three components (kinetic, ther- 
mal and magnetic, respectively) of the total energy in each shell: 

Ek (F, p. Ax) : = F(F - 1 )pc'^ Ax 
ETiT,p,p,Ax) := l(pe + p)T^ - p]Ax ^^^^ 

EM{T,p,cr,Ax) := |f^ - ^jpcrc^Ax 

where p, p, e, a and F are the fluid rest-mass density, thermal pres- 
sure, specific internal energy density, magnetization parameter and 
the Lorentz factor. Ax is the width of the shell. Then, the total dy- 
namical efficiency is defined as the ratio of the sum of the magnetic 
and thermal energies in the post-shock state to the total initial en- 
ergy (MAIO). 

Recently, NKTll have computed the radiative efficiency of 



magnetized internal shocks using a slightly different definition. The 
definition of NKTl I has the advantage of resulting into a radiative 
efficiency, which is Lorentz-invariant, while the one of MAIO is 
not. Therefore, the dynamic efficiency can be computed in the lab 
frame and then used in the rest frame of the contact discontinuity 
to be directly compared with the radiative efficiency, as was done 
in Sec.|5] However, NKTll approach does not shed any light on 
the problem of obtaining a true energetic efficiency in terms of the 
initial conditions of the shells. This is because their definition does 
not consider the efficiency of conversion of the initial total energy 
of the shells to radiation but, instead, the efficiency of conversion 
of the enthalpy per particle after the shell collision into radiation 
enthalpy. 

With the aim of introducing a Lorentz invariant, energy-based 
definition of the radiative efficiency in terms of a quantity akin to 
the dynamic efficiency of MAIO, we consider the following expres- 
sions for the kinetic, thermal and magnetic energies 

Ek(T,p,Ax) := FV^Ax 
Et{T,p,p,Ax) := r^(ps + p)Ax . (A2) 

£m(F,P, cr. Ax) := F^pcrc^Ax 

With these definitions, the ratio of energies shown in Eq. [3T] 
becomes Lorentz invariant, and so is the total dynamic efficiency. 
Furthermore, when calculated in the LAB-frame, the exact values 
of the dynamical efficiency computed using any of the two sets 
of definitions (Egs. lAll or IA2t differs very little, as it is shown in 
Fig.lAl] 



APPENDIX B: COMPUTATIONAL REQUIREMENTS FOR 
REALISTIC CALCULATIONS 

In this paper we have adopted a cylindrical geometry and ig- 
nored adiabatic losses of the non-thermal electrons, as well as 
the high-latitude emission. However, we have included both of 
these effe cts in previous works studying afterglows of gamma- 
ray bursts jMimica et al .|2010| :Im imica & Gianniosl20T l|). In those 
calculations, perform ed using the radiative transfer code SPEV 
jMimica etani2009l) . we have only taken into account the syn- 
chrotron and EIC emission processes. Even so, the calculation of 
a realistic multi-wavelength light curve of a single model lasts 
anywhere between few hours to few days on a supercomputing 
cluster with several hundreds of computing cores, depending on 
the number of wavelengths at which the emission is computed. If 
the synchrotron- self absorption (SSA) is included, than the paral- 
lel scalability of the method is reduced and the ca lculations can 
become prohibitively expensive (see footnote 3 of iMimica et al.l 
I2OIOI) . 

In this paper we have a similar problem due to the fact that we 
include the SSC process, which, as SSA, is non-local and difficult 
to parallelize in a realistic conical geometry. The use of cylindrical 
geometry and other simplifying assumptions has enabled us to nev- 
ertheless compute light curves of 900 models for 96 frequencies at 
120 observer times. This required 200 thousand computing hours 
on the MareNostrum computer of the Barcelona Supercomputing 
Center. Based on our previous experience, we estimate an order of 
magnitude more resources are needed if the adiabatic cooling of 
the electrons is included, and another order of magnitude if an ac- 
curate radiative transfer method such as SPEV is used instead of 
approximate method described in Section |4] We note that such a 
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Figure Al. Left panel: contours of the total dynamic efficiency for the model parameters from Table 1 computed using definition of MAIO (see Eq. lAll 
Contours indicate the efficiency in percent and their levels are 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 and 13. In the region of the parameter space above the dashed 
line there is no forward shock. Filled contours show the magnetic efficiency. Right panel: same as left panel, but using Eq. lA2l inspired by NKTl 1. 



calculation requires resources in the range of 10 million computing 
hours on a computer such as MareNostrum. 
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